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The B-mode polarization lensing signal is a useful probe of the neutrino mass and to a lesser extent 
the dark energy equation of state as the signal depends on the integrated mass power spectrum 
between us and the last scattering surface. This lensing B-mode signal, however, is non-Gaussian 
and the resulting non-Gaussian covariance to the power spectrum cannot be ignored as correlations 
between B-mode bins are at a level of 0.1. For temperature and B-mode polarization power spectra, 
\£) , the non-Gasussian covariance is not significant, where we find correlations at the 10 -5 level even 

for adjacent bins. The resulting degradation on neutrino mass and dark energy equation of state 
is about a factor of 2 to 3 when compared to the case where statistics are simply considered to be 
£S) , Gaussian. We also discuss parameter uncertainties achievable in upcoming experiments and show 

that at a given angular resolution for polarization observations, increasing the sensitivity beyond a 
certain noise value does not lead to an improved measurement of the neutrino mass and dark energy 
equation of state with B-mode power spectrum. For Planck, the resulting constraints on the sum 
, of the neutrino masses is o"E m „ ~ 0.2 eV and on the dark energy equation of state parameter we 

find, a w ~ 0.5. 

PACS numbers: 98.80.Es,95.85.Nv,98.35.Ce,98.70.Vc 
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OS . I- INTRODUCTION 

The applications of cosmic microwave background (CMB) anisotropy measurements are well known its ability to 
constrain most, or certain combinations of, parameters that define the currently favorable cold dark matter cosmologies 
with a cosmological constant is well demonstrated with anisotropy data from Wilkinson Microwave Anisotropy Probe 
""^5 ' • Furthermore the advent of high sensitivity CMB polarization experiments with increasing sensitivity [j| suggests 
that we will soon detect the small amplitude B-mode polarization signal. While at degree scales one expects a unique 
.B-mode polarization signal due to primordial gravitational waves Q , at arcminute angular scales the dominant signal 
will be related to cosmic shear conversion of -E-modes to B-modes by the large-scale structure during the photon 
propagation from the last scattering surface to the observer today 

This weak lensing of cosmic microwave background (CMB) polarization by intervening mass fluctuations is now well 
studied in the literature 0, Q , with a significant effort spent on improving the accuracy of analytical and numerical 
calculations (see, recent review in Ref. [8J). As discussed in recent literature the lensing B-mode signal carries 
important cosmological information on the neutrino mass and possibly the dark energy, such as its equation of state 
0, as the lensing signal depends on the integrated mass power spectrum between us and the last scattering surface, 
weighted by the lensing kernel. The dark energy dependence involves the angular diameter distance projections while 
the effects related to a non-zero neutrino mass come from suppression of small scale power below the free-streaming 
scale. 

Since the CMB lensing effect is inherently a non-linear process, the lensing corrections to CMB temperature and 
polarization are expected to be highly non-Gaussian. This non-Gaussianity at the four-point and higher levels are 
exploited when reconstructing the integrated mass field via a lensing analysis of CMB temperature and polarization 
[lOj. The four point correlations are of special interest since they also quantify the sample variance and covariance 
of two point correlation or power spectrum measurements [Tlj . A discussion of lensing covariance of the temperature 
anisotropy power spectrum is available in Ref. ^2] • I n the case of CMB polarization, the existence of a large sample 
variance for B-modes of polarization is already known [l3l | , though the effect on cosmological parameter measurements 
is yet to be quantified. Various estimates on parameter measurements in the literature ignore the effect of non- 
Gaussianities and could have overestimated the use of CMB -B-modes to tightly constrain parameters such as a 
neutrino mass or the dark energy equation of state. To properly understand the extent to which future polarization 
measurements can constrain these parameters, a proper understanding of non-Gaussian covariance is needed. 

Here, we discuss the temperature and polarization covariances due to gravitational lensing. Initial calculations on 
this topic are available in Refs. while detailed calculations on the CMB lensing trispectra are in Ref. [l^. Here, 

we focus mainly on the covariance and calculate them under the exact all-sky formulation; for flat-sky expressions 
of the trispectrum, we refer the reader to Ref. [T^j. We extend those calculations and also discuss the impact on 
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cosmological parameter estimates. This paper is organized as follows: In SjnJ we introduce the basic ingredients for 
the present calculation and present covariances of temperature and polarization spectra. We discuss our results in 
£ 11111 and conclude with a summary in flIVI 



II. CALCULATIONAL METHOD 



The lensing of the CMB is a remapping of temperature and polarization anisotropics by gravitational angular 
deflections during the propagation. Since lensing leads to a redistribution of photons, the resulting effect appears only 
at second order [8|. In weak gravitational lensing, the deflection angle on the sky is given by the angular gradient of 
the lensing potential, S(h) = V</>(n), which is itself a projection of the gravitational potential <f>: 

0(m) = -2 r\r dA{r y r) M r,rhr), (1) 
J Q dA(r)dA{ro) 

where r(z) is the comoving distance along the line of sight, r$ is the comoving distance to the surface of last scattering, 
and dA(r) is the angular diameter distance. Taking the multipole moments, the power spectrum of lensing potentials 
is now given through 

{^Im&lm) = 5ll'8rn m 'Cf (2) 



as 



Cf = \j k 2 dkP(k)l}™(k)lf»(k), (3) 

where 

pro 

I l r(k) = \ drW Xen {k,r)h{kr) , 
Jo 

W^k,r) = -^n(^) 2 F(r) d fy r \ , (4) 
V k ) dA{r)dA(ro) 

where F(r) = G(r)/a(r) and G(r) is the growth factor, which describes the growth of large-scale density perturbations. 
In our calculations we will generate Cf^ based on a non- linear description of the matter power spectrum P(k). In the 
next three subsections we briefly outline the power spectrum covariances under gravitational lensing for temperature 
and polarization E- and S-modes. In the numerical calculations described later, we take a fiducial flat-ACDM 
cosmological model with fl b = 0.0418, Q, m = 0.24, h = 0.73, t = 0.092, n s = 0.958, A(k a = 0.05 Mpc" 1 ) = 2.3 x 10~ 9 , 
m v = 0.05 eV, and w = — 1. This model is consistent with recent measurements from WMAP 0- 



A. Temperature anisotropy covariance 

The trispectrum for the unlensed temperature can be written in terms of the multipole moments of the temperature 
0i m as 15] 

+ Ch c lA- l ) mi+m2 Kh mi5 "X m ' 3 ■ (5) 

It is straight forward to derive the following expression for the multipole moment of lensed 9 field as a perturbative 
equation related to the deflection angle 0: 

9~lm = 0l m + ^ fcmi^m2^™ imi + 2 <t>hrm #Z 2 m 2 <t>l 3 m 3 ^U^C^^ ' ( 6 ) 

l 1 m 1 l2m 2 l\in\l 2 ra 2 l 3 m 3 

where the mode-coupling integrals between the temperature field and the deflection field, I™^ 1 ™ 12 and "/i™^;™ 2 ™ 3 , 
are defined in [l5l[Tr3|. 
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As for the covariance of the temperature anisotropy powerspectrum, we write 

C ° Vee = 2h + 1 2h + 1 S ^hmAmJhmJt 2m2 ) ~ Cf^ = O + V + (Q + tt)* llIa (7) 

mim,2 

where the individual terms are 

P = (2h + lK2i 2 + l)E C M^Fui,Fm 
4 \^ n<t>r e r e (i? \ 2 

- t(2h + l) ^ ° tiC ^ fL{L + 1){2L + 1] ' (8) 



Q 



and the last two terms, which are related to the Gaussian variance, can be written in terms of the lensed temperature 
anisotropy power spectrum as 



2 

Q + K 



^nr, (9) 



where 



cf = [i-a 2 +o^+E^%TT c f 2 

h 



- >(i, + i) + i,(l, + i)-.(i + i)V (a + 1)ph ; 1)(afa + 1) (i J do) 

We note that Eqs. HlOfl are readily derivable when considering the lensing effect on the temperature anistropy spectrum 
as in Ref. @. 



B. i5-mode Polarization Covariance 



Similar to the case with temperature, the trispectrum for an unlensed £7- field can be written in terms of the 
multipole moments of the £?-mode Ei m : 

(E hmi E hm2 E hm3 E hm4 ) = ^^(-lJ^^M^^H^C^-lJ^^^Cr 4 

+ QfQ?(-i) mi+m2 cr m4 cr m3 - ( n ) 

To complete the calculation, besides the trispectrum of the unlensed E- field in Eq. Ijllfl. we also require the expression 
for the trispectrum of the lensing potentials. Under the Gaussian hypothesis for the primordial i?-modes and ignoring 
non-Gaussian corrections to the cf> field, the lensing trispectra is given by 

+C? x Cf 2 (-l) mi+m2 ^™7 m4 <5 ; X m3 • (12) 

For simplicity, we assume that there is no primordial B field such as due to a gravitational wave background and 
find the following expression for the lensed E'-field: 

E lm = E lm + \ ]T ^ imiJ B i2m2+2 /;rj 2 lim2 (i + (-i) Wi+i2 ) 

+\ e ^ imi ^ 2m ^?3 m3 +2^r 2m3 ( i + (-i) i+h+h+i *) , (is) 

Iimi(2m2l3m3 
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FIG. 1: The correlation matrix [Eq. 1251 1 for temerature (left), B-mode (middle), and B-mode (right) power spectra between 
different / values. The color axis is on a log scale and each scale is different for each panel. As is clear from this figure, the 
off diagonal correlation is weak for both 9 and B-mode power spectra, but is more than 0.1 for most entries for the B-mode 
power spectrum. This clearly shows that the non-Gaussianities are most pronounced for the B-mode signal and will impact the 
information extraction from the angular power spectrum of B-modes than under the Gaussian variance alone. The B-mode 
covariance shown in the left panel agrees with Figure 5 of Ref. 14] . 



where the expressions for the mode coupling integrals +2 / / ™'^ im2 and +2Ju™ 2 \™ :trn3 are described in Refs. p^LITr|. 
As for the covariance of _E-mode powerspectrum, we write 

Covee = 2l l - 12l l +l E (E hmi E: imi E l2m . 2 E; 2m2 ) - Cfdf = H + 1+(J + JC)S llh (14) 

^ m 17712 

where 
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1)(2Z 2 H 
2 


-1) 


{2h + 


1)(2Z 2 H 
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(2Zi + 


1)2 E 
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C$1 


(h(l 


1 + 1)- 


4) 



3 = : 1 , 9 >;c^^ , (i+(-i)' i+l+ '')(2^ 1 l/0 2 



E C t( C hf L ( L + l)(2i + 1) . (15) 



2 7 r(2Z i 



The last two terms can be written in terms of the lensed power spectrum of -E-mode anisotropies as 

2 



^ = — (C^, (16) 



where 



8tt 

h 



2-Rl,i a = ![,,(,, + 1) + m + 1) - ,(, + I)] / 2 ' + '" 2 '- + ™»> + " %) . (17) 



Note that Cf is the power spectrum of the lensed i?-modes. 
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FIG. 2: Here we show the cumulative signal-to-noise ratio for a detection of the power spectrum [Eq. 1311 1 for temperature 
(left), _E-mode (middle), and B-mode (right) polarization power spectra. The solid line is the case with a Gaussian covariance 
whereas the dashed line is with a non-Gaussian covariance. We can see that for the case of the temperature and E-mode 
polarization there is little difference between the Gaussian and non-Gaussian covariance, but for the B-mode polarization there 
is a difference of a factor of ~ 10 at large I values. 



C. B-mode Polarization Covariance 



The calculation related to B-mode power spectrum polarization is similar to the case of the i?-modes except that we 
assume that the B-mode polarization is generated solely the lensing of the i?-mode polarization. Based on previous 
work (c.f., Ref. Q), we write the multipole moments of the lensed -B-modes as 



iBi m = \ £ iimi ^ 2 „ l2+2 /™J 2 lim2 (l-(-l) Wl+i2 ) 

l\m\l 2 m 2 

+ \ E ^m 1 S I> m 3 ^3 Tn ,+ a TO, mam, ( 1 -(- 1 )' +,1+,a+,S )- ( 18 ) 



4 

(imii2m2(3m3 

Here, we will only calculate the £?-mode trispectrum with terms involving Cf since we will make the assumption that 
corrections to B-modes from the bsipectrum and higher-order non-Gaussianities of the lensing <fi field are subdominant. 
Thus, using the first term of the expansion, we write 

(B hmi Bi 2m2 B hm3 B hmi ) =— Y, E E E 

LiMil^m^ L 2 M 2 l 2 m' 2 L 3 M 3 l' 3 m 3 L^M^l'^m'^ 
{(f>L 1 M l 4'L 2 M 2 (t)L 3 M 3 4'L 4 M i ) + 2ll imi L! 

(19) 

where 

+2li mhmi i 2m2 = 2F lh i 2 (-l) m ^ m ^ ^ J (20) 

The covariance of the B-mode angular power spectrum can be now defined as 

Covbb = 2/ +121 + 1 ^ {Bi xmi Bl imi Bi^ m2 Bl 2m2 ) - Cf x Cf 2 ■ (21) 

After some straightforward but tedious algebra, we obtain 

Covbb = A + B + C + 8 hh V, (22) 
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FIG. 3: The derivatives of the temperature (8), _B-mode, and B-mode power spectra with respect to the sum of the neutrino 
masses (oc Q,„h 2 , top panel) and the dark energy equation of state, w (bottom panel). It is clear that in the case of the sum 
of the neutrino masses the addition of the B-mode polarization greatly increases sensitivity. In both cases we find that large I 
information also increases sensitivity. We note that the derivative of the temperature power spectrum with respect to neutrino 
mass agrees with that shown in Fig. 3 of Ref. [T^ . 



where the terms are given by 



A = 



4(2Zi + l)(2Z 2 + l) 



E 



(Cj) 
2L + 



2 / h+L 



i E ^(i-(-i) h+i+i ')( 2 ^^) 2 



. l'=\l x -L\ 



h+L 

E C^l-{-lt +L+l ')( 2 F l2Lll 

A> = \l 2 -L\ 



r, N E 
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4(2li + l)(2I 2 + l)^ 
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FIG. 4: The expected error on the sum of the neutrino masses (top three panels) and the dark energy equation of state, 
w (bottom three panels) as a function of experimental noise for three different values of the beam width, #fwhm- The solid 
line considers Gaussian covariance with just temperature information, the dotted line considers non-Gaussian covariance with 
just temperature information, the dashed line considers Gaussian covariance with both temperature and polarization (E- and 
B-mode) , and the dot-dashed line considers non-Gaussian covariance with both temperature and polarization. It is clear that as 
the beam width is decreased the estimated error on the sum of the neutrino masses and w is increasingly overly optimistic when 
just the Gaussian covariance is used in the Fisher matrix calculation. We choose 5 bins uniformly spacing between I = 5 and 
/ = 100, while we choose 13 bins logarithmic uniformly spacing between I — 100 and I = 2000. This choice of bins are sparser 
compared to |2(i|. From the expressions of covariance matrix [Eas. fft . 11411 . 1221 1 . we know the Gaussian parts are diagonal and 
therefore the larger the bin is, the more important the non-Gaussian effect is. So the non-Gaussian effects in our bandpower 
statistics are more obvious that those in (2y|. 



2 N <t> h+Li h+Li h+h r £ ^ - 

C = 16(2/ 1 + l)(2/ 2 + l) E E E E C L G t G i C l [2F hLlVl 2F hL2V22 F hLlV22 F hL2Vi ] | £ * * ^ 



where 

c? = lE c t (j wr^ cg(i - (-l)^). (24) 

Unlike the calculation for the covariances of the lensed temperature and polarization i?-mode, the numerical cal- 
culation related to covariance of the B-modes is complicated due to the term C, which involves a Wigner-6j symbol. 
These symbols can be generated using the recursion relation outlines in Appendix of Ref . , though we found that 
such recursions are subject to numerical instabilities when one of the I values is largely different from the others and 
the I values are large. In these cases, we found that values accurate to better than a ten percent of the exact result can 
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FIG. 5: The error ellipses from our Fisher matrix calculation. We have varied eight parameters, and show the error ellipses 
for each parameter with Em„. The dot-dashed ellipse is the expected error from Planck with just a Gaussian covariance, the 
solid ellipse is same but with a non-Gaussian covariance. The short-dashed ellipse is for an experiment with the same beam 
width as Planck (6 ~ 5') but with decreased noise (1 fjK^/sec as opposed to 25 /^K v / sec) with a gaussian covariance and the 
long-dashed ellipse is the same but with a non-Gaussian covariance. 

be obtained through semi-classical formulae 0|- I n am7 case, we found that C is no more than 1% of A, B, and these 
terms are in turn no more than 10% of T>. The same situation happens to those expressions in flat sky approach 20J. 



III. RESULTS AND DISCUSSION 



We begin our discussion on the parameter uncertainties in the presence of non-Gaussian covariance by first estab- 
lishing that one cannot ignore them for the B-mode power sptectrum. In Figure 1 we show the correlation matrix, 
which is defined as 

Cov XY (i,j) (25) 



y/Cov X Y(i,i)Cov X Y(j, j) 
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This correlation normalizes the diagonal to unity and displays the off diagonal terms as a value between and 1. 
This facilitates an easy comparison on the importance of non-Gaussianities between temperature, E-, and i?-modes 
of polarization. As shown in Figure 1, the off diagonal entries of temperature and E- modes are roughly at the level of 
10~ 5 suggesting that non-Gaussian covariance is not a concern for these observations out to multipolcs of 2000 [l^ . 
while for B-modes the correlations are at the level above 0.1 and are significant. 

Below when we calculate the signal to noise ratio and Fisher matrices, we use the bandpowers as observables with 
logrithmic bins in the multipole space. Our bandpower estimator for two quantities of X- and Y-fields involving 
temperature and polarization maps is 



A 2 



1 «« 1 / 

ly yi 

C*i ^ ^— ' 47T 

I— In m=—l 



Xl m Yr 



Iw 



(26) 



where oti — la — In is an overall normalization factor given by the bin width. The angular power spectra are 

Af = (Af) = ^^ + l)ZCf^, (27) 



while the full covariance matrix is 



with the Gaussian part 



((A 2 -A 2 )(A 2 -A 2 )HS«% + ^, 



S,;,; — 



(28) 



(47r) 2 a 2 ^ 
' i — tii 

A 



MB 



i(* + l)0FWH M /81n2 



and the non-Gaussian part is 



S N = - 



5^(2Ii + l)(2l 2 + l)hh{Cov^ Efi ) . 



(29) 



(30) 



hh 



To further quantify the importance of non-Gaussianities for B-modes, in Figure 2, we plot the cumulative signal- 
to-noise ratio for the detection of the power spectra as a function of the bandpowers. These are calculated as 



(SN) 



Cf i r Cov^(A„A J )C 



XY 



(31) 



by ignoring the instrumental noise contribution to the covariance. As shown, there is no difference in the signal- 
to-noise ratio for the temperature and -E-mode power spectra measurement due to non-Gaussian covariances, while 
there is a sharp reduction in the cumulative signal-to-noise ratio for a detection of the £>-modes. This reduction is 
significant and can be explained through the effective reduction in the number of independent modes at each multipole 
from which clustering measurements can be made. In the case of Gaussian statistics, at each multipole I, there are 
21 + 1 modes to make the power spectrum measurements. In the case of non-Gaussian statistics with a covariance, 
this number is reduced further by the correlations between different modes. If N is the number of independent modes 
available under Gaussian statistics, a simple calculation shows that the effective number of modes are reduced by 
[1 + (N — l)r 2 ] when the modes are correlated by an equally distributed correlation coefficient r among all modes. 
With N = 21 + 1 and substituting a typical correlation coefficient r of 0.15, we find that the cumulative signal-to-noise 
ratio should be reduced by a factor of 7 to 8 when compared to the case where only Gaussian statistics are assumed. 
This is consistent with the signal-to-noise ratio estimates shown in Figure 2 based on an exact calculation using the 
full covariance matrix that suggests a slightly larger reduction due to the fact that some of the modes are more 
strongly correlated than the assumed average value. 

To calculate the overall impact on cosmological parameter measurements using temperature and polarization spec- 
tra, we make use of the Fisher information matrix given for tow parameters [i and v as 



p — 



E 

X = B,E. 



e ij 



i 

XX 



d(Af) 2 
— a — ■ 



(32) 
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where the summation is over all bins. While this is the full Fisher information matrix, we will divide our results to 
with and without non-Gaussian covariance as well as to information on parameters present within temperature, and 
E- and B-modes of polarization. 

Since B-modes have been generally described as a probe of neutrino mass and the dark energy equation of state, 
in Figure 3, we show dCi/dm v and dCi/dw to show the extent to which information on these two quantities are 
present in the spectra. It is clear that B-modes are a strong probe of neutrino mass given that the sensitivity of 
temperature and B-modes are smaller compared to the fractional difference in the B-modes. Furthermore, B-modes 
also have some senitivity to the dark energy equation of state, but fractionally, this sensitivity is smaller compared 
to the information related to the neutrino mass. 

In Figure 4, we summarize parameter constraints on these two parameters as a function of the instrumental noise 
for different values of resolution with and without non-Gaussian covariance. While for low resolution experiments 
the difference between Gaussian and non-Gaussian extraction is marginal, non-Gaussianities become more important 
for high resolution experiments where one probes B-modes down to large multipoles. In this case, the parameters 
extraction is degraded by up to a factor of more than 2.5 for both the neutrino mass and the dark energy equation 
of state. We have not attempted to calculate the parameter errors for experiments with resolution better than 5 
arcminutes. This is due to the fact that such experiments will probe multipoles higher than 2000 and we are concerned 
that we do not have a full description of the non-Gaussian covariance at such small scales due to uncertainties in the 
description of the matter power spectrum at non-linear scales. As described in Ref. @, the CMB lensing calculation 
must account for non-linearities and their importance only become significant for small angular scale anisotropy 
experiments. Furthermore, we also do not think any of the upcoming B-mode polarization experiments with high 
sensitivity, which will be either space-based or balloon-borne, will have large apertures to probe multipoles above 
2000. 

The value of 2000 where we stop our calculations is also consistent with Planck. Since Planck HFI experiment will 
have a total focal plane polarization noise of about 25 /zK -^sec, based on Figure 4, we find that it will constrain 
the neutrino mass to be below 0.22 eV and the dark energy equation of state will be determined to an accuracy 
of 0.5. Note that the combination of Planck noise and resolution is such that one does not find a large difference 
between Gaussian and non-Gaussian statistics, but on the othcrhand, experiments that improve the polarization noise 
well beyond Planck must account for non-Gaussian noise properly. In future, there are plans for a Inflation Probe 
or a CMBpol mission that will make high sensitive observations in search for a gravitational wave background. If 
such an experiment reach an effective noise level of 1 [iK ^sec and has the same resolution as Planck, the combined 
polarization observations can constrain the neutrino mass to be about 0.18 while the dark energy equation of state 
will be known to an accuracy of 0.44. This is well above the suggested constraint from Gaussian noise level. This 
suggests that while high sensitive B-mode measurements are desirable for studies involving the gravitational wave 
background, they are unlikely to be helpful for increasingly better constraints on the cosmological parameters. 

The non-Gaussianities in the B-modes, while providing information on gravitational lensing, limits accurate pa- 
rameter estimates from the power spectrum alone. This is contrary to some of the suggestions in the literarture that 
have indicated high precision of measurements on parameters such as the neutrino mass and the dark energy equation 
of state with CMB B-mode powerspectrum by ignoring issues related to non-Gaussian correlations. Furthermore, 
while atmospheric oscillations suggest a mass-squared difference of Am^ ~ 10~ 3 for two of the neutrino species, it is 
unlikely that one will be able to distinguish between mass hierarchies with CMB polarization observations alone if 
one of the two masses related to the atmospheric oscillation result is close to zero (c.f. 19] ). This is discouraging, 
but understanding the information present in CMB polarization beyond powerspectra, such as direct measurements 
of non-Gaussianities themselves, could potentially allow an improvement. 

From Figure 4 we see that as we decrease A p the measurement errors on the parameters asymptote to a constant 
value. We can understand this in the following way. As we see from Eq. I|29|) . the noise blows up exponentially at large 
I and therefore sets an effective cutoff ^o. Only the bandpowers which are smaller than Iq contribute to parameter 
estimates. Therefore, if we decrease A p , we increase the number of bandpowers we can observe and hence obtain 
better sensitivity with negligible instrumental noise for I < ^o- Therefore, the curves in Figure 4 become flatter as 
we decrease A p . The same situation applies to Figure 5. Figure 4 also shows that as we decrease the beam width, 
$fwhm, we see the Gaussian covariance becomes more significant. This is a result of the fact that the Gaussian 
covariance grows in significance with increasing I. 

In Figure 5, to highlight the impact on cosmological parameters beyond the neutrino mass and dark energy equation 
of state, we also show constraints from the Fisher matrix calculation. We show error ellipses calculated with and 
without the non-Gaussian lensing covariance for two different experiments: Planck, with 6*fwhm = 5' and A p = 
25 /iKy^sec and 'super-Planck' with $fwhm = 5' and A p = 1 /iK^/sec. This comparison shows that while parameters 
such as m„ and w are affected, parameters such as r, fl m h 2 are not affected by non-Gaussian information. This is due 
to the fact that the cosmological information on these parameters come from temperature and B-modes rather than 
B-modes. This highlights the fact that the issues discussed here are primarily a concern for the B-mode measurements 
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and extraction of parameters, especially the parameters that have been recognized to be mostly constrained by the 
B-mode measurements, and not for temperature and E'-modes. 



IV. SUMMARY 

The i?-mode polarization lensing signal is a useful probe of certain cosmological parameters such as the neutrino 
mass and the dark energy equation of state as the signal depends on the integrated mass power spectrum between us 
and the last scattering surface. This lensing £>-mode signal, however, is non-Gaussian and the resulting non-Gaussian 
covariance to the powerspectrum cannot be ignored when compared to the case of temperature and polarization E- 
mode anisotropy covariances. The resulting degradation on neutrino mass and dark energy equation of state is about 
a factor of 2 when compared to the case where statistics are simply considered to be Gaussian. We discuss parameter 
uncertainties achievable in upcoming experiments. 
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